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f# • iti 


Itiaaticn C3«@9) gi^aa i 


i3»m} 


0 • 0 • al lli« liiili 

• 1 § a* wia nog# isMi sMiwittasx aaionup 

Btffareniiatdi^ •gpaatiait C3«3f ) uilii faaanag In « and itatcn 
atpalian C3«3)» wa gall 


f a • («, . to4 ^ » (lb • K«> d* 

IQlt 
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C3#30l 


MXm^ 


fsilig <3*2 )t 
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unMiP nil %r 
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(3.33) 


ttoo IliiMENiSL ^usitaRSf laopw 

•fOttkian |a* 13 l ii iwsi aai ti^Liii llii C 3 «i 9 > t« 

C3«3t> «• iMiiilr 




i 




^tt ©*4 




0 « © 
0*1 


■Yj « 0 


Cili tuft - \ 

Clwlfl 

C«fe %)!« «c{i« of tilt 
liioiml ^rnm^mf lagror i 


§iiie« mm listoroit l4o« la aodarato oiwaitaffo# 
wita tiii ppaetiaal viisf folafet m mil ^meoaferolMi mm arnlymt 
far falocitjr tHofMO. IwwBtor 1 j«w* fwr aotoato ©arfmaror 
’I?Smi Hinwil tNiuatoiy lifar C3«B1> fm 

•mrfuittNi mail isoiaiiary mmMttmm CS*S^> i» I^MKrafor* to 
imlrod* HIM vftlaoa of ^ irnmiMamiaMil ilroia iwmiitm f - 
oi mmlmm oalito* ^ rj mm mmi«d ia to otoii • ioimioa 
for tliia oottotioa* 

too ymlmtlm lafwr owtian CS*1il 

aaOor tuo lM»tiiteif oonOitioiii BpmiMmi lioa olroilir 

Hoiii oolvoa fir Cff t*» f**lli!r Itofloiia osamri C1»i)t m 
Umm mtmA m aoto iio|im«tla:itm aofcisod for tto ooliimM oo 



f 4l tWbolir ioiMiovr Saoror Bfo^ioat 

’ftpftMMi 10 a imrnm ooOor aoamdtaoMr 

mfMmsfW MMtmmMtA mmUmm it it rnfftmill to 



ixstd trlaJl aotf mtrm 

m»Wm4. ffitmtlar l&o fwt « ftited m3m of paxrm^m 

te# Milottfd f f***<0) th« tsolp of wofiiwi loioif* 

1he 

l®4i« to 8toi*t the integrotioti tt the ^»r 

li«i% eiMpj^ed mii^e of f*ff** i^re eoioeroi homiiefar 
ooiiilitiixie* Besei oa' the reet^tiog ^iwe «t the outer IMlt* 
nm Vetoes of f lO) wii f |0> mm 9m0m nM iitegn^loo 
mmtMrnim Iho wmomm woe eoittifiitiKl teitil tho 

eoeorac^ ^of f * at the outer hoonteip wm hotter thei 
8 X 10 me Hhit of t** hetl«r tho» j t x 10 • 

’Roege * Mto etof hf eter oethM of intiirelieo wee tieei* 
ih liPL. the liioio ooimi% tt wee iteeteeeif to etort 
the ihiegfotiiii hr oelo^tnc waloee of mikmm Hsitiiil. eonii* 
tiQii% heiei on jphreloet tiilotioii m oi^ieni&ao» t» fM 




aiiorooiid^l® mtmt in wiliiiiwt* 

.In to imid diotroi^tiitios tiiMn*<int in 

tiitt sSiGm eiitli»(l»» C3«12) iMioii oolvnd tiiii^ 

toelani^ vbio'h Mo i»«nki Mplainod ia 

Bin i^*liii«»r nation CS«^2I io of Hm Hcmt 


f**** nfCTi.ftfifiV**) «0 a. 3 < 

•Hi non iMi nrittiii S» ttio fbvw •• timm In 



»id© 0f tli« Amm tli# 

£ioii«^aaaii«i}«mi« Suffis n itimttcMiw 

Sli« imiiiidiiT w«iclitif%i9 tli^ CS*37) 


«r«t 

f «0) « 0 

f *(0) - 0 

f'C^a} ■ »"■ ' {"■ ' » at tiM atiga of tlia vatoeltgr toonaafy 

C1«4'^^F ir®lo@ity 


at« two eimdltioBo g** wi g*** ar« n*^ tetioMii m% rj m % 

I«no« Wmm will too two oc^lom^t^fr i^latloiso mA mm purt^* 
emXm intoftmi m lotpliiliioi to tppooitx«€« . 

fHi tnitiol ooi^ttoiig timM io mim&AMAmm it 
« 0 fpr tuo furot oiwiplliiietixr oolutloii It of oipKlt®tt 
Witt mMHmmmAmmm tpmg to »«po ovo to to 
toliiii tot 


f (0>ot 
g*^ CO) m 0 

r” (0) ■ f 

ir’**(o) • 0 


(3.M) 


fto imiioaL imSitiQiii g«r sflooiid oonitiottiiitoff mhMm M 

#g ofMlIiMi CI«3gl ollli iioii«4tii«oitii«^ ^mm ppiol to «o»o 


mm to to IuOmb oos 


t (0) • 0 
t* ( 0 ) • 0 
«*' CO) ■ 9 


(S.)» 



rang®® tmm « 0 at to® surfan® t® “Vj • 7*1 at tli® 

9 &$& of tn® ]^®r» Sli® irs^Ui® of 7| at to® oif® of 

to^ bomitetT lairor oii^ to tor uiiito is siioifiototlr 
iafgo* It tel® tisto lOisid toat 7| m 7*1 is too most siaitslil® 
oteloo* toito if o® ^>00® tor nito toto® is a® 

ai^ooitol® oliaiii® la to® toltoiio* 

to® Mipi.it® tol^ioa of «®aitito 0#371 sito 
lioisitofr osiliit i W Bi® (3*8t) 'tot f3*8!fJ is ttoi iiaoto soifci®®** 





sf' C^#*! • li| <K JSI^ ♦ Cfji H iSSt ♦ SI’S* 

f**|?#i> ■ ®| » «!*• ♦ ©n * a** ♦ «•* 0.>f5l 


SeMng C 3 «M^> an© O^h ^ ^ 

©I 

©i^yoi Wm inl.ii«» &t ©I •»© ©n tteui e^taio«i aai 
tit© «© 34 itliiiia fit tt© ^ ^i«t 9 ©tviiwliv©© ff*» 

»% «s*» it**# a**# a**» a***t a***© a***© ««© tai«*g 

©f f© f*© f**© f*** •*• ©I ©UrfMili gl•tl©Cl» 

f *>1 • ©» 1 © ©•ft « # « t© 7 ] • f»®)* 

«s©©ft©ltt ifiMtiM* Um ©©3itii©ii Cf© f*» f**© f*’*) 
©Inm© Ik© Ifftiniitf «i ^*wfe t«*it t, « i i a# ©oaiSliMnilttiV' iMift 

i!MhllP©ii» tiij^ <iiiajiifetai>©^ t© ©t| © ©©©©©fit IPI©* Hi© ©MMtiPiJfiajL ©pl^tt 

Myjj# afiA^ ||l | | ^ aiMi OiMMUPH^ 

1^ mwimA mUsMm ©t ©©Hi it©ti«ii Cv|« Dbf© ^•t© %% 

* • « tt'Vj « f*PI* a Hmi itfltrpM# t«|ii©©ii Wm mmmmM 

•^IgliMi PMI Hm iMir fplHliM «t mw alMlM p4fit it MVt 



CB) ro'JncUST E^tiam 

tilt tlians^ boim^ikXT lasrar C3*33l 1® % 

Xlufiir tiOKagiiiaoiit action af mmM afdwr irith 
aaniillMsIStM)* I3i# iiimX mtimd af %li« Itft f<»iii% 

llna^ Baiif^Lwqr $f*fliaiaR ii m9d» ttm .protadisra rm§ m 
falterat 

ftpatioit C3«B3) i« riwrittaii mi 

0** ♦ ,.lfcJC,J.! « & (l«33l 

aii4A I 


in# tiiilial. ami tiawitoy aanaitiaiia aval 







« 0 11,% rnsettm® m% tim «%• of t&i« 

lat^or^ %i3@ rmwm imt uMioli Ims air«i»f33f liooii uliitoa* 

?ii« s^l^idii of «<pK%liiKi C3*33l i« ^ 

Ui}««r «cMi%%ii&ti<iis of Mt m4. Z& is 0mm 

0 «i 0 « lai ♦ mt C3«MII 

^mm 0 it • 0iiis%«us% Mxm t» to l>« msii^ Mi# 

#at«r iMwiOajqr ^loAitioii C 0 « 1 «l It 

X# ofdir t# olufft fmt OOtaiiJtiiig %!» aoliMoiit Mi# 
firtt ioitiilxifi &• ##tiin«0 «# %# mm Iftyil i# 

C# • 0 t 0*« 01 Hi #11 tiM 
'^l • fttit 

f iliii Ml# 93mm toiiikion (t « 0«. &*m 01^ Mi# 

%if# lift n OUtaliiMI idMl IMIittl MiiiiiM# 



ilwffi Wsm iTiilne of %iit l> • Otiiig 

of Sllf ZM o nm eol^tifin it 

.iOiiorioA mioos tuts nm ioiiiilaii uro oMixurod 
i#illi IteiO of %hm tolotto ot ^mU oi^iMost iioint* "Qm 

pfoeooii if ropoitwl ifi^ nm mhittm m mmirvjl ooitit^ioii 'till 
4iislr®d aoctiaraosr if ;iolii@trot« 

Xi is to to iioto4 tliat tlio tisii*4liia«:isi«(fNA strcMWi 
f«etioii to oqoitloii Cl* 33) if ttm tli® mluttm 
of tte oolooitr boif^ofy Issfor o^tiso C3*IS^ «o4 ttm HiorwA 
U^w i» solwa for 10*®*) far f rani® of 

i*r»oatl muhmrm 

maferilmtlQB .aoiiif^liQQi 

4 t to® Odj^cit of 9«wr»tio% til# otiftir iireti i« mto 
ttist if f** o 0* flii® prorMo® m oO^itionil initial oon^tion 

»'* m mO 

111* vwstuamrlaal wloel^ tomiMT ( 3 « 3 r> i* 

no® solfoi in tiio mmm olroa^ 4»«®riMl wuSm mo flolloiflni 
initial mA tomiiiT o®o4iti.oiioi 
f C0)o0 

t* (0) > 0 
f” CO) *0 

r* >-JU. «» •«■• of «b» taOBilMT )ar«r 

14*V1 : 


fti9 mlmtion tfUl Ktmwim of a mh^im 

tnitiM. o^mMitods 
f 0} 
f* CO) • 0 
i** imrnQ 
f***CO)»t 

and a pitftie^lAr adtuHoct xit^ initial. 
t CO) • 0 
f* CO) • 0 

f** <0) > 0 
f’" (0) « 0 

fHo final anlnticHt fm strain fnmstimi f# elitalsiai 
«t sainimtiiii lua bane issad to find tnmrattira diatrUmtiQn 
at sopiiAitioit for tlia tliaroal hmm&mr l«f«r CS«3S) 

tHa CB«3f)« naM mmirsnal tn^Micpt i» 

rapaatai in Umi idicnm mntioiiod inl«f ratitoa* 



mmm m 


EssoTO Mm mmmmtm 


Wafm® for dtsemssicm m mtmiX 

mA l«i«p«ra.tiiFo dtstflbutioii* fior otarroA rnMtmm^ II is 
tt«irii3»I# lo Ills mQwmw of mnsriml. Isilsiiiw 

wm Ea'so follcMsS* Hiii iimt Iomi dioo# for tim fliMOlal. tm«« 
of s flfti f 3.»|« for Mtm Wm '©orfmlaro 4 « 0 • '&• mlmtW 
iiilrilii^im rssoili lisiro tuMjklti .iii M&lm CX) 
plolloiS nfsliisl' r| $M. ftgtir# Cf ) for loaPiotm iloroiiimft m 
omaoiitod In soellon 3t^« XI U mm tln^ tmfilii iloro* 
tifln giroo losi rsmlls tinioii oonpuPtn mitf i^xx tno 
wvmllsilil# Slitliu solnlioii for liis rsloeilf for n 

nnl fXnlo* m mmm fmm ftpiro IHo iroloiii^ 

•r«Si«ii| f** ol lOo won nijiiiiod If ilio mrnimd roiis 0«I3I 
ns sfsiiisl • dtllE 19 ^ iln^nin 



vsjMitviP msirijppiiiOB Bwii wSiB onisiiinioa iw n 


f«i^ <if snrvnliivo wsnslsr | 4 « « 0*01 lo s 0*01 I §m 


mwimm vnlnsi of r| C 'Vj •» § to 7| • I* Hio foti^io 

99W MM UPl# Mm pwmIMI Hi |#lii 


mlmiXf for ©'usirattu?® Iiai atlm 

hmn AoMi miffing ftm tii« point Cii/Q| « 1» ^ ■ 0)*. 

Iipotion 1m9 itood to pIM potonti^ 

'diitn^i^^ioisi* '%# iSn® ?9 Join® ®XI tli® pe^Ms tli® 

potmtiaX dlrtritutioii Mno® mort tli# mtmitj 4i0* 

tfitutiois Mu#® iit tl}« tiou^dajnr r®gi^ for ttfformi 
oiirtmtiroi* fliii iin® ro^«a«ot® tte liftit cMT tli# 

tliieJffio®® tm" musA oiafirod mmtmp* It i® i^on tti^l tM» 
tomdaiT TlMg'm ttiio*^®® ineroM®® 4® turvmtiar® inoroM®® 
froK nogfttiir® to poittii^® Tmiwi « i^oolt 
toot ottaSnot tgr pf««ioa» w»rM«r®^\ 

fjitl® (XI*4l iiidioat#® ro'loisi^ OirtrlliittiiMi^ at 
iMipHmtitit for rartoiia <wr®ttxr®« at difr®r««it r\ « 

'Hi® ronaits aro .pliitt®il in ifU Sia« wplmlW it vm 

m% til® Mrfiir® aiii tit® fotintiai roXoeitjr at It® ^ 

®«li» tto® iMOwiiaigr lartr at i® ^ oat® for Hi® v®]®®iir 
iiatrslitttitii im Wm tcmitrr Xs^r®r rflgt<m mam M J^Lgiir® (#)• 


t fli® r®i^ of ourvtiitr® pyriMittr 4 li®® tNiin ototo® 
ti«Hiiiii 4 • • t® 4 « ♦ ®»0| l»«®«it®« it I4®« in Hi® 

rang® ®f frartittil. iiit«r®®t« f®r •xaig»l®t flgiir® C4> it®®® 
tiii^ tm 4 • ® i«@it Hi® iTiiiii® ®f ^ at lisi® mtMP ®®i® ®f ^ 




Mlm tfm iraaisfH • 0 is Hi tliis e4»» 

'Xliis mmns th^ ttm 3tm^ irtrssi at thm iioiiit of 

soiKirationi is wm&* It nan bo i»m flptros 

CSI iniiS C?lt timt ttio mignitiMo of iroloeit^ at angf 7| tmt m 
partl«a.lar ewv%%wt9 it isoro btforo soorirmtioii m it ibouM 
b<i« 

Xt is to bo mt&^ tem Fiftiroi CS) C7)t that 
^ tlio offaot of ounratoro on wXooilf aLstfibtatiofi i» mm 
prononteod at iraltDis of ^ 2 that it mw fbem ttio titrf«ieo» 

Hi® .profU® got® «o4i.fi®S for iafi® Vj 4oo to tb® 

ottiar boi»4of|r owiiitioo* 

ta fifofo Ci)» tb® offset of owiratiif® i® ulotta® 

Ob tb® tabgiotiol tolooitf ifilimt f** at tii® siixfiO® (rj m Oh 
fitit Hmir® that tb® sbMr itr®®® for oomtoa sorfa®®® is ^l®«i 
tlisit tit® fiat plat® mill® vliil® that for ooiiear® siprfaooa i® 
iroator ttaii Hi® fl^ plat® mlo®* Xt mmmt that tli® ocmni 
rarfi®® is son ssbadtin*® to a® aivirs® isroMbOf® sribiiobt 
(2 C '^/ §3 Hiwi Wm iBammm smtam9 M olb®fiiofd% 

for a oopfwK «iif®®®« Hi® Hioir tHr«®® iNiiiii losst ^ mmm^ 
tiob Mila laM® piM® ®®riiar if tbof® i® mi aorori® prsiaiiNi 

* |(!n 

'™(w ms" Up 

fimraiar® iistnimtio® iia® b«®B cmliiilmi im 
mmimm mx»®® ®f omptaiitm wrmMtor C a ® » dsCii » a • 0 
mi a « # ibil' ) «*' MStmmk v] ^ Hi® rtisAi® bsMi inmi 
t» Mi®® |xix» m«at ^ •»& xf«al m wmm% 






ilititliir. is ists suolitous isssisss tnsrs is iis» assd 
sf iusssii^ ills ioiiisi ssn^itisii sisli %$m^. MSm 
Hi# ^msMl^sjr #<i»diMiG» #1 Hi# ouisr #%• #f Wm 
itfjr is ful^^r mtlsf i«d is Hi# ursssuft 

m» rssnits wm in irsxy^ i##«l nftsuMit mIHi 

Hist# tfP Hsssi^^^* 

f* wrfms^ #f mirvstiirtt m Ussi trsssfar %t#s iartwaili* 
iiisi tm sisrfasss wim f«itiis ^ ewt 9 m%'m^ 
t# sfw# toon #f instep# tn# slaiiiiiii#ii piinit 
liliit Hi# isstM ^#r stellar ##liill#iis» Hi# li«i» Hmot* 
ter »#i# i# »#r# ter • #«i@«f# tmtmm mm ter Hi# 
#####• mxnmm mm rmm^ is stellar te hsHsp# «te* 
SliiiiteA W PWSiM rsribirs fsr s^ fsteteolte 
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TABLE n I A 

TEMPERATURE DISTRIBUTION FOR VARIOUS CURVATURES 

AT PRANDTL M0=1.0 
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temperature distribution for various CURVATU 

AT PRANDTL N0=1.0 
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TABLE IV 


tlmplrature distribution For various curvatures at separation 

AT PRANDTL NO=1.0 
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TABLE IV A 

temperature distribution for various Curvatures AT separation 
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appendix d 


COMPUTER PROGRAMME 


$JOB MEG147 * TImE008»PAGES 010 *NAmE K N SRiVASTAVA 
siBJOB Map 

SIBFTC MAIN NODECK 

dimension Y(91» 5»4)*Z2(91» 5»4),Z1{?1» 5 »4 ) ♦ X ( 9 1 ) » Y 1 ( 91 ) » 

1 Y2(91 ) »Y3( 91 ) »Y4(91 ) »Y5{91 ) *Z3(91» 5*4) 

2 .P(91» 5,2)»ZZ2(91* 5»2)*ZZ1(91, 5 , 2 ) »P1 { 91 ) »P2 { 91 ) » P3 ( 91 ) 

N=NUMBER OF STATION POINTS 
ETA infinity TAKEN EQUAL TO 7.1 
equivalence {Z1>ZZ1 ),(Z2,ZZ2) 

N»81 
A=-0»06 
DO 36 111=1,7 
A=A+0.02 
PRINT11,A 

11 FORMAT (2x,F8, 2) 

C ASSUMED SOLUTION 

DO 6n0 1 = 1 ,N 
Y( 1,1,1 )=*0.0 
Y(I,1,2)=0.0 
y( I ,1,3 )=n.o 
YC I ,1,4) =0.0 
600 CONTINUE 

17o K=0 

K = K+1 

C COM'^L SOLUTIONS Zl Z2 

yl(l)=0.0 

Y2( 1 )=0.0 
y3(1)=0.0 

yA(1)«1,0 I 

Y5(1)»0.0 I 

DO loo I=2,N 

IF( I.GT.ll ) GO TO 222 

H=O.Ol 

GO TO 333 

222 H=0.1 , 

333 call RUNGA(Y1( I-l ) ,Y2 ( I-l ) ,Y3 ( I-D ,Y4 ( I-l ) ,Y5 ( I-l ) ,0*0 ,Y ( I ,K 

1 Y( I ,K,2 ) ,Y { I ,K ,3 ) ,Y ( I »K,4 ) ,Z1 ,X , I ,A ,K+1 ,H ) 

y1(I)=X(I) ! 

Y2( I )=Z1( I ,K+1,1) ! 

Y3( I )=Z1( I ,K+1,2) 

YA(I )=Z1{I,K+1,3) ' j 

y 5( I )=Zl( I ,!<. + l,A) I 

loo continue 

[ Yl(l)=0.0 

I y2(1)“0*0 

j Y3(1)=0#0 i 

1 YA(1)«0*0 ■ . 

i Y5(1)=1.0 ' 



DO 2 Of) 1=2, N 
IF(I,GT.11) GO TO 444 
H=O.Ol 
GO TO 555 

444 H ” 0 • 1 

5 55 call RUNGA(Y1( I-D ,Y2 { I-l ) ,Y3 ( I-l ) ,Y4 { I-l ) ,Y5 ( I-l ) » 0. 0 , Y ( I ♦ K , 1 ) » 

1 Y( I ♦K,2) ,Y( I ,K,3) ,Y( I ,K,4) ,Z2 ,X,I ,A,K+1 ,H) 
y 1{I )=x(n 
y2(I )=Z2( I,K+1,1) 

Y3( I )=Z2( I,K+1,2) 

YA( I )=Z2( I ,K+1,3) 

Y5(I )=Z2( I,K+1,4) 

200 CONTINUE 
c particular integral 
Y 5(1 )=0.0 
DO 4o0 1=2, N 
iFd.GT.ll) GO TO 666 
H=0.0l 
CO TO 777 

6 66 H = 0*1 

777 call RUN6A{Y1( I-l ) »Y2 ( l-l 1 ,Y3 < l-l ) »Y4( I-D , Y5 ( I- 1 ) » 1 -O * Y ( I ,K ♦ 1 ) * 
1 Y( I .K»2 ) »Y( I »K»3) I »K»4) ♦Z3,X,I ,A,K+1,H) 

Yl ( I ) =X { I) 

Y2( I )=Z3( I ,K+1,1) 

Y3{I )*Z3( I,K+1,2) 

YA(I )=Z3( r,K+l,3) 

Y5(I )=Z3( I ,K+1,4) 

A-00 continue 

Y(1»K,1)»0.0 

c USE OF third And fourth Boundary conditions 

B»1.0/( 1*0+A*7,1) 

BB«-A/( ( 1.0+A*7,1 )**2) 

ZlP»Zl(N*K+l,2) 

Z2F*z2{Nti<+l,2) 

Z3P=Z3(N,K+1 ,2 ) 

Z1DP=Z1{N,K+1,3) 

Z2DP=z2(N,K+1,3) 

Z3DP»Z3(N,K+1,3) 

c 

c calculation of constants 

c . , : 

cl = ( B*Z2DP-BB*Z2P-Z3P«Z2DP+Z3DP*Z2P ) / ( Z1P*Z2DP-Z 1DP*Z2P ) 
C2={B-Z3P-C1»Z1P) /Z2P 
DO 3oo 1=2, N 
DO 3 oo J“1*A 

Y( I,K+l*J)=Cl*Zl( ItK+l*J)+C2*Z2(I .K+lfJ)+Z3(I»K+l,J) 

3 00 CONTINUE 

C TESTING THE ACCURACY 

Y(1»K+1,1)=0.0 
Y(1»K+1,2)«0.0 
Y(i»r+i»3)»ci * ; 

Y ( 1 tK+l ,4 ) »C2 ^ 


■I 



K1=K+1 

PRINT 211 ,(Y(M,K1,2),M=1,N) 

211 FORMAT (12X,10(F12. 8)) 

DELTA=0.0000001 

1=2 

180 IF(ABS(Y(I,K+ 1 * 2 )-Y<I,K* 2 ) UlT. DELTA) Go TO 99 

DO 77 I=1,N 
YM,K,1)=Y{I ,K+1,1) 

Y( I.K,2)=Y(I ,K+1,2) 

Y{ I ,K,3)=YM *K+1,3) 

Y(I,K,4)=Y(I ,K+1,4) 

77 CONTINUE 

GO TO 170 
99 1=1+1 

IFCI.EQ.N) go to 190 
GO TO 180 
19o X(3)=0.0 

DO 35 1 = 1, N ,10 
X( r)=Yl(I) 

J = K+1 

Y( 1»J»3)=C1 
Y{1,J,4)=C2 

PRINT26»X(I) ,Y(I,J,1) ,Y(I,J.2) fY(I»J,3) ♦Y(I»J,4) 

26 F0RMAT(2X,5F15,8 ) 

35 -CONTINUE 

C 

c 

c energy equation 

c 

PR=0.0 

do 701 1 1=1 #2 
CALL FLUN (20000) 

PR=PP+0.5 
PRINT 111, PR 

111 FORMAT (/20X,F15.8) 

DO 6ol I=1»N 
P( I,1,1)=0*0 
601 P( I ,1,2 )=0.0 
171 k: = 0 

X = K+1 
Pl(l)»0.0 
P2(l)»o*0 
P3(l)=l,o 
DO lol I=2,N 
IF( I.GT.ll) GO TO 881 
H=0.0l 
GO TO 991 
881 H=0.1 

991 call RUNG(P1{I-1) ,P2 ( I-l ) ♦ P3 ( I -1 ) ,0 . 0 .P ( I »K , 2 ) , PR . Y ( I , J * D , 
1 ZZ1,X*I ,K+1 ,A,H ) 

Pl( I )=X( I ) 

P2(I)=ZZ1(I*K+1»1) 

■p:3T':i:):*zz ' - 



P3(l)=0.0 
. DO 2ol r=2,N 
IFd.GT.ll) GO TO 888 
H=0.0l 
GO TO 999 
888 H=0.1 

999 call RUNG(P1 Mr-l) ,P2( I-l) , P 3 { I -1 ) , 1 • C »P ( I »K » 2 ) » PR » Y ( I , J ♦ 1 ) ♦ 
1 ZZ2,X» I ,K+1,A,H ) 

Pi ( I )=Xt I) 

P2(l )=ZZ2{ r,K+l,l) 

2(1 P3( I )=ZZ2( I,K+1,2) 

P{1,K,1)=0.0 

C=(1.0-ZZ2(N»K+1 »1 ) )/ZZl(N»K+l»l) 

DO 3ol I»2,N 

p( i ,K+i»i)=c*zzi( r»K+i,i)+zz2( 

301 P(I,K+1.2)=C-»-ZZ1(I,K+1,2)+ZZ2(I,K+1»2) 

P( 1,K+1*1 )=0.0 
P(1,K+1,2)=C 
1=2 

181 IF(ABS(P( I ,K+l«l)-Pn ,K»1 ) )*LT.DElTA) GO TO 133 
DO 78 1=1, N 
P(I,K,1)=P(I ,K+1,1) 

P{ I ,K,2)-P( I ,K+1,2) 

78 CONTINUE 

GO TO 171 
133 1=1+1 

iFd.EQ.N) GO TO 191 
GO TO 181 

191 no 38 1=1, N, 10 

JJ=K+1 ' 

PRINT 27,X( I ) ,P{ I ,JJ,1) »P( I »JJ»2) 

27 FORMAT (2X,15F8. 4) 

38 CONTINUE 

70l continue 

36 CONTINUE 

STOP 
END 



C RUnGE-KUTTA integration subroutine For velocity profile 

c 

sibftc 

SUBROUTINE RUNGA { R1 ,R2 , R3 , R4 ,R5 , DEL ,SF ,S1 ,S2 ,S3 , Z ,X »N , V, I 1 *H ) 

dimension A{A) ,B(A) »C( A) ,Y(6) ,Q ( 10 ) »DY ( 6 ) ,2 { 91 , 5,A),X(91: 

A{1)=0.5 I 

A(2)=l,-(0.5)**,5 

A{3)=l,+ (0«5)-»-*,5 

A(A)=1, 0/6.0 

BID =2.0 

B(2)=l.o 

B(3)=l.n 

P{4)=2.0 

C(l)=0*5 
C{2)=li-SQRT( .5) 

C(3)=1.+SQRT{.5) 

C( A)=0.5 

Yd) =R1 d 

Y(2)=R2 
Y(3)=R3 
Y(4)=R4 
Y(5)=R5 
DO 25 1 = 1,, 5 
Q( I )=n.n 
25 continue 

50 DO lo J=1 

DYd )=1,0 I 

5=3 .O+V* ( Y( 1 )+H/4.0) 

DY(2)=Y(3) i 

DY(3)=Y(A) I 

DY(A)=Y( 5) I 

DY(5 )=-“( ( 2.#V*S3 ) /S+V#-»-3*S1/S**3-V«-V*S2/S#»2+ S1*S2/ ( 2 .*S )+SF*: 
12.4S)+V*SF*S2/ (2.*S**2 ) +V*S1*S1/ ( 2 .*S»*3 ) -V*V*S1* ( SF-Sl* { Y ( 1 ) +! 
2{ 2.*S**3 ) )*DEL-( Y{2)-SF*DEL)*(S3/(2.*S)+V»S2/(2.*S**2)-V*V* S| 

3 .*S»*3)) -(Y(3)-S1*DEL)*(V»*3/S»*3+S2/(2.«S)+ V«-S1/{S)**3 

4 -{SF#V*#2-2.#V*V*S1»(Y(1 )+H) )/(2.*S«-»3> ) 

5 -(Y(4)-S2*DEl)*(-V*V/S**2+S1/(2.*S}+V*SF/(2.*S**2) ) 

6 -.(Y(5)-S3*DEl)*(2.*V/S+SF/{2.»S) ) 

•fDO 2o 1=1,5 ■ ; 

temp=a( j)*{Dy( I )-B( j)»q{ n ) 

Y( I )=Yd )+H*TEMP 

Q( I ) =Q( I )+3.0*TEMP*C( J)»DY( I ) 

20 CONTINUE 

10 CONTINUE 

X(N)=Yd) 

Z(N,I1,1)=Y(2) 

Z(N,I1,2)= Y(3) 

Z(N,I1,3)=Y(A) 

Z{N,I1,4)=Y(5) 

RETURN 

END 



c rUnGE-kUTTA integration subroutine for temperature profile 
c 

sibftc eng 

subroutine rung (R1*R2»R3»DEL,YY *PR*YP »Z,X»N» II ,V,H) 

dimension A(4) ,B(4) »C{ A) ,Y(6) ,0 ( 10 ) ^DY ( 6 ) » Z ( 91 , 5»4)»XI 

A( 1 )=0.5 

A{2)=1,-(0.5 )**.5 
A( 3) =l.+{0.5 )**.5 
A(4)=l, 0/6,0 
B( 1 ) =2.0 
B{2)=1.o 
B(3)=l.n 
B(4)=2.o 
C( 1)=0.5 
C( 2)=1.-SQRT( .5) 

G{3)=1,+SQRT(.5) S 

C(A)=o,5 

C H=0*1 . i 

Y( 1 ) =R1 ! 

Y(2)=R2 
Y(3)=R3 

DO 25 1=1,3 f 

25 Q(I )=0.n i 

50 DO lo J=1,A . I 

S“1.0+V»(Y( l)+H/4.0) i 

DY(1 )=1.0 ’ 

DY(2)=Y(3) 

DY(3)=-(PR-*-YP*YY/(2.0*S) )*DEL 
1 -{Y(3)-YY*DEL)*PR*YP/(2,0»S) 

DO 2n 1=1,3 

TEMP=A( J)*(DY( ! )-B( J)*Q( I ) ) 

Y{ I) =Y( I 1+H»TEMP 

Q< I )=Q{ I )+3.0*TEMP-C( J)»DY( I ) 

20 rONTiNUE 

10 CONTINUE 

X(NT = Ytl> : ■ 

■ Z(N,I1,11*y:{2) - S'K'-.." V: 

Z(N,I1,2)=Y(3) 

RETURN I 

• " :END ^ 

sentry ■ i 
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